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ABSTRACT 

Aims. One of the major uncertainties in close binary evolution is the efficiency of mass transfer /?: the fraction of transferred mass 
that is accreted by a secondary star. We attempt to constrain the mass-transfer efficiency for short-period massive binaries undergoing 
case A mass transfer. 

Methods. We present a grid of about 20,000 detailed binary evolution tracks with primary masses 3.5-35 M G , orbital periods 1-5 days 
at a metallicity Z = 0.004, assuming both conservative and non-conservative mass transfer. We perform a systematic comparison, 
using least-squares fitting, of the computed models with a sample of 50 double-lined eclipsing binaries in the Small Magellanic 
Cloud, for which fundamental stellar parameters have been determined. About 60% of the systems are currently undergoing slow 
mass transfer. 

Results. In general we find good agreement between our models and the observed detached systems. However, for many of the 
semi-detached systems the observed temperature ratio is more extreme than our models predict. For the 17 semi-detached systems 
that we are able to match, we find a large spread in the best fitting mass-transfer efficiency; no single value of /? can explain all 
systems. We find a hint that initially wider systems tend to fit better to less conservative models. We show the need for more accurate 
temperature determinations and we find that determinations of surface abundances of nitrogen and carbon can potentially constrain 
the mass-transfer efficiency further. 

Key words, binaries: close - Stars: evolution - binaries: eclipsing - binaries: spectroscopic - Stars: fundamental parameters - Stars: 
mass loss - Magellanic Clouds 
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^ 1- Introduction 



Evolutionary calculations of massive close binaries have 

been conducted with v arious ste llar evolution codes 

since the 1960s, e.g. Paczvnski (f]966l) . 



(I 1 9691). iDe Greve & De Loord 



Kippenhahnl 
11994 . 



"(11 992b. IPols 

De Loore & Vanbeyerenl (1 19941) . IWellstein et al.l (1200 lb and 
Nelson & EggletorJ (12001 1) . One of the major uncertainties 
in these calculations is the efficiency of mass transfer: what 
fraction /3 of the transferred mass is actually accreted by the 
secondary star? Conservative evolution, i.e. the case where mass 
and angular momentum of the binary system are conserved, 
corresponds to (3 — 1. Non-conservative evolution has been 
considered in most calculations by assuming a constant f3 < 1 
while the amount of angular momentum loss is described using 
a second parameter. 

This question is relevant for many current astrophysical 
problems as it affects, for example, the for mation rate of progen - 
itor models of long-soft gamma -ray burst (Petrovic et al.ll2.Q05b , 
of double neutron star binaries dDewi et all 20061 2007 ). which 
are thought to be the progenitor s of short-duration gam ma-ray 
bursts, and Type la supernovae (lYoon & Langerll2005l) . which 
are used as standard candles to measure universal expansion. 

Two effects of mass transfer on the accreting star are likely 
to resul t in ma s s loss from the system: expansion and spin up. 
iBensonl dl970l) . lUlrich & Burgerl d 19761) and others have shown 
that when mass transfer occurs on a timescale on the order of 



the thermal timescale of the donor star, the less massive accret- 
ing star, with a longer thermal timescale, is driven out of ther- 
mal equilibrium and expands. In systems with initial mass ratios 
very different from unity, this leads to a contact or common- 
envelope configuration. Significant mass an d angular momen- 
tum l oss from the system is then expected dFlannery & Ulrich 
1 19771). The second effect, spin up, was pointed out by iPacket 
(Il98ll) : after gaining only a few percent of its original mass 
through disk accretion, enough angular momentum is transferred 
to spin up an isolated accreting star to critical rotation. In prin- 
ciple this can lead t o sign ificant mass loss from the system (e.g. 
lLanger et alJl2003l 120041) . depending on how efficient tidal ef- 
fects can keep the accreting star rotating in synchrony with the 
orbit. 

No consensus has been reached on this topic. Many relevant 
processes are still not yet well understood. Hydrodynamical sim- 
ulations are the most promising approach, but are still far too 
time-consuming to study the dependence of the mass transfer 
efficiency on the binary parameters. At the moment the most 
fruitful way to address the efficiency problem is to parametrize 
the process of mass transfer and to use observations to calibrate 
the parameters. The most stringent tests come from double-lined 
eclipsing binaries for which the stellar parameters can be deter- 
mined with accuracies of a few percent. Semi-detached systems 
undergoing their first phase of mass transfer provide the best test 
objects, to avoid uncertainties introduced by an unknown mass 
loss history. 
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Various studies in which theory is compared to observations 
indeed indicate evi dence for non-conservative mass transfer. 
iRefsdal et alj ( 1 19741) showed that it is very likely that the semi- 
deta ched b i nary A S Eri is the result of non-conservative evolu- 
tion. [Sarna|([l993|) studied the semidetached system /3 Per (Algol) 
and claimed that it has lost about 15 percent of its initial total 
mass and 30 percent of its initial total angu lar momentum. Also 
for P Lyr moderate mas s loss w as inferred (iDe Greve & LinneiH 
1 1994b . iFTgueiredo et alj d 1994ft compared Galactic OB binaries, 
of which 8 are semi-detached, to stellar evolution models. They 
estimate that these systems have lost between 30% and 60% 
of the transferred mass. Many Wolf-Rayet binaries with O-type 
companions require a highly non-conservative firs t mass-transfer 
phase to explain their orbital periods and masses dPetrovic et alj 
120051) . IVan Rensbergen et al.l d2006l) compared a grid of calcu- 
lated binary models statistically to observed orbital periods and 
mass ratios of Algols. They need to assume a significant amount 
of mass loss to obtain agreement between models and observa- 
tions. 

Other studies, however, show tha t mass transfer is fair l y con- 
servative, at least for some systems. iNelson & Eggleton! d200 lb 
compared observations of semi-detached binaries with inter- 
mediate mass, in which both components have spectra in the 
range GO to Bl, to models assuming conservative evolution and 
they found an acceptable agreement overall. However, they did 
not compare the observed systems to non-conservative mod- 
els. Almost conservative evolution during the first phase of 
mass transfer is also needed to explain the formation of sev- 
eral types of evolved bi naries, including so me massive X-ray 
binaries, e.g. Wray 977 dKaper et alj 1 19951) . and eccentric bi- 
nary systems consisting of a white dwarf and a neutron star 
dPortegies Zwart & Yungelsodll999tlTauris & Senneisll2000h . 

In view of the observational evidence for non-conservative 
mass transfer, in many calculations of binary evolution a sin- 
gle constant mass transfer effici ency has been assumed for 
all systems often J 3 = 0.5 (e.g. IDe Greve & De Loorelll992t 
IChen & H an 2002). However, neither the above-mentioned ob- 
servational comparisons nor the theoretical considerations out- 
lined earlier provide a basis for such an assumption - in fact, 
they clearly indicate that /3 is not a constant but probably de- 
pends on the masses and orbital properties of a binary. No clear 
picture of this dependence has emerged yet, partly hampered by 
the fact that much of the evidence is based on incidental studies 
or heterogeneous observational samples. 

In this paper we attempt to shed more light on this ques- 
tion by means of a systematic comparison between binary evo- 
lution models and a ho mogeneous sample of binaries undergoing 
mass transfer. Recen tly lHarries et alj |20()3. hereafter |H03|) and 
iHilditch et al.l d2005l hereafter |H05|) presented the fundamental 
stellar parameters for a sample of 50 double-lined eclipsing bina- 
ries in the Small Magellanic Cloud (SMC). It is the largest single 
set of fundamental parameters determined for high mass stars 
in any galaxy. More than 50% of the systems are semi-detached 
and believed to be currently undergoing the slow phase of case A 
mass transfer (mass transfer in a system consisting of two main 
sequence stars). As case A mass transfer is the first phase of mass 
transfer after the stars left the zero-age main sequence, there are 
no uncertainties induced by previous mass transfer phases. The 
detached systems in the sample enable us to test our models in 
the pre-mass transfer phase. 

No suitable set of case A binary models are available at this 
m oment. Large grids have b een calculated before, for example 
bv lNelson & Eggletonld200ll) and lWellstein et al.l d200ll) , but not 
at the metallicity of the SMC, i.e. Z = 0.004. With this work we 



present a large grid of detailed case A binary evolution models 
at the metallicity of the SMC with different assumptions for the 
mass transfer efficiency. We address the problem of mass transfer 
efficiency by fitting binary evolution tracks to ea ch i ndivi dual 
system in the observed sample presented by |H03l and lH05l We 
then investigate if correlations can be found between the best 
fitting mass transfer efficiency parameter and the initial binary 
parameters. 

2. Stellar evolution code 

To calculate detailed binary evolution tracks, we used the 
STARS stellar ev olution code, a varian t of the code origi- 
nally developed b y lEggletonl dl97llll972l) . An important update 
dPols et alJI 19951 and references t herein) was the impro vement 
of the original equation of state (Eg gleton et al.l [19731) by in- 
clusion of pressure ionization and Coulomb interactions, OPAL 
opacity tables and recent nuclear reaction and neutrino loss rates. 
A recent addition to the code, which we have used in the calcu- 
lations presented here, is the so-called TWIN mode in which the 
structure and composition equations for both stars in a binary 
are solved simultaneously with equations for the spin and orbital 
angular momentu m, the orbital e ccentricity and the mass flux 
between the stars dEggletonl200 6f1. 

Convective mixing is modeled by a diffusion equation for 
each of the composition variables. A mixing length ratio l/H p = 
2. is assumed. Convecti ve overshooting is taken into account as 
in ISchroder et al.l d 1997b with a free parameter i5 ov =0.12 cali- 
brated against accurate stellar data from non-interacting binaries 
dSchroder et al.lll997hlPols et al.lll997l) . 

In close binary systems tidal interaction tends to circular- 
ize the orbit and to synchronize the orbital period with the 
spin periods of the stars. Spin-orbit interaction by tides is 
treated according to the equilibrium tide theory dHutl Il98lt 
Eggleton et al.lll998b . We enforce synchronized rotation in all 
our models by decreasing the tidal friction timescale by a 
factor of 10~ 4 with respect to the default value given by 
lEggleton & Kiseleva-Eggletonl d2002l) . The effect of rotation on 
stellar structure is taken into account as a reduction of the effec- 
tive gravity assuming rigid rotation. 

The prescription for mass transfer implemented in the TWIN 
version of the code allows us to treat both semi-detached and 
contact binaries. Mass transfer by Roche-lobe overflow is mod- 
elled as a function of the potential difference S between the stel- 
lar surface and the Roche-lobe surface, for each star. As long as 
neither star exceeds its Roche lobe (tf> s < 0) the mass-transfer 
rate is zero. If one of the stars overfills its Roche lobe, the mass 
transfer rate is calculated by solving an additional differential 
equation for the mass flux at each mesh point outside the Roche 
surface potential, simultaneously with the other structure equa- 
tions, 

^.-Cx^, (1) 
dm r 

where m is the mass coordinate and r the radius. The mass trans- 
fer rate is then given by the integral of eq. (Q~|i over all mesh 
points outside the Roche surface potential. By choosing a suffi- 
ciently large value for the free parameter C, i.e. 10 4 , we make 
sure that the radius stays close to the Roche lobe radius in semi- 
detached binaries and the mass transfer rate is self-regulating. In 



1 A description of the current version of the code can be obtained on 
request from ppeOigpp .ucllnl . org. 
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the case of stable mass transfer the mass loss rate does not de- 
pend on this function, it is set by the evolutionary expansion rate 
and the change in orbital separation. The exact form of eq. <£TJ 
is therefore not important. The main rationale for using it is that 
we can treat the case where both stars overfill their Roche lobes 
in a similar way. The direction of mass flow now depends on the 
difference in surface potentials <p s \ - <p S 2 between the stars, 



dMi 
dm 



= -C x sign(f sl - (j> s2 ) 



V 2 l<^i - <M 



(2) 



with the equivalent equation for star 2 having the reverse sign. 
Hence mass flows from star 1 to star 2 if (p s \ > </> S 2 and vice 
versa. 

We note that the physics and evolution of contact binaries 
is a complex and currently unsolved problem in stellar evolu- 
tion. Massive contact binaries are quite common, and although 
they are less well studied than their low-mass equivalents, the 
WUMa systems, both classes appear to share the property of 
having components with nearly equal effective temperatures. 
This requires a mechanism for heat transfer between the stars. 
A simple physical model for heat transfer in low-mass contact 
binaries with convective envelopes was developed and included 
in the TWIN version of the code by Yakut & Eggleton (2005). 
However, the physics of heat transfer in binaries with radiative 
envelopes is probably quite different. For want of a quantitative 
model, we ignored the possibility of heat transfer between the 
stars in the calculations presented in this paper. 

Mass and angular momentum loss during mass transfer is 
modeled by treating the mass transfer efficiency as a free param- 
eter, 



P = -MaIM d , 



(3) 



where A stands for accreting star and D for the mass donor. j3 is 
assumed to be constant throughout the evolution. We assume that 
the mass lost from the system takes away the specific angular 
momentum of the orbit of the accreting star, 



h = — — a a>, 



(4) 



where M tot represents the total mass of the system, a the sep- 
aration between the two stars and a> the orbital angular veloc- 
ity. Reverse mass transfer, i.e., mass transfer in an evolved sys- 
tem from the initially less massive component to the initially 
more massive component, is assumed to be conservative. This 
assumption does not affect the comparison with the observed 
sample as all stars are on the main sequence and reverse mass 
transfer occurs only in evolved systems. 

The rate of stellar wind mass loss decreases with metallicity 
dVink et al.ll200lHMokiem eta l. 2006). For comparison with the 
observed systems our main interest is the core hydrogen burning 
phase during which the amount of wind mass loss is small com- 
pared to later stages of evolution and compared to the amount of 
mass loss during Roche lobe overflow. Therefore we assume no 
mass loss in the form of a stellar wind. 

We calculated our models at a metallicity of Z = 0.004. For 
the hydrogen and helium abundance we assume 



X = 0.76-3.0Z, 
Y = 0.24 + 2.0Z. 



(5) 



The abundances of the heavier elements are assumed to scale 
to solar and meteoric a bundance ratios as determined by 
lAnders & Grev esse (1989) so that they are consistent with the 



opacity tables (R ogers & Iglesiaslll992|) even though the actual 
abundan ces in the SMC deviate to some extent from solar pro- 
portions dRussell & Dopitall 19921) . 

3. Binary evolution tracks 

We have calculated four data cubes of binary evolution trackfl 
with different initial primary masses Mi,;, mass ratios q, = 
M[j/M2j and separations or orbital periods P,. The subscript ; is 
used to indicate that these values are the initial values, i.e. at the 
moment both stars start burning hydrogen on the main sequence. 
These parameters are spaced at approximately equal logarithmic 
intervals, 



log Mi, = 0.55, 0.60, 1.55, 
log 4, = 0.025, 0.050, 0.350, 
log(P,/P ZAM s) = 0.05, 0.1, 0.75, 



(6) 



with Pzams an approximation of the orbital period at which the 
initially more massive component would fill its Roche lobe on 
the zero-age main sequence for a system with equal masses, 



ZAMS 



0.19Mi,, + 0.47Mf: 
1 + 1.18Mf ( . 



2.33 



(7) 



One data cube of binary tracks is calculated assuming con- 
servative mass transfer (J3 = 1). Three additional data cubes 
are calculated assuming non-conservative mass transfer (J3 = 
0.75,0.5 and 0.25). Each binary evolution track is terminated 
when one of the following conditions is met: 

- 4000 time steps are taken, 

- The stellar radius exceeds the Roche lobe radius by more 
than 10%, which happens if deep contact is established, a 
situation our code is not designed for. 

- The mass loss rate becomes larger than 3 x 10 2 M/tkh, where 
M is the mass of the initially more massive star and tkh its 
Kelvin Helmholtz timescale. These systems are expected to 
come into contact. 

- The code fails to converge. This occurs for example during 
advanced burning stages and in models for which the mass 
loss rate becomes very high. These models are expected to 
come into contact shortly after. 



3. 1 . Subtypes of case A evolution 

Our conservative grid of models is very similar to the survey 
of case A binary evoluti on at solar metallicity published by 
iNelson & Eggletonld200ll) . In line with their definitions we iden- 
tify different subtypes of case A binary evolution. As an exam- 
ple the subtypes are plotted in Figure[T]for two different primary 
masses and for different assumptions about the mass transfer ef- 
ficiency p. 

- Case B: In our widest systems the initially most massive star 
will not fill its Roche lobe before depleting hydrogen in the 
core. A system is classified as case B if it fills it Roche lobe 
after leaving the main sequence, which we define as the mo- 
ment the central hydrogen abundance drops below 0.001. 

- Case AR Contact during rapid mass transfer: In systems 
with large mass ratios, the mass-accreting star reacts to the 
mass exchange by expanding rapidly as it is driven out of 



2 Our models are available on-line from the following website 
http://www.astro.uu.nl/data/stars/ 
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Fig. 1. Classification of different types of binary evolution. In each diagram the sub-type is given as a function of initial orbital 
period and mass ratio for a fixed primary mass and fixed mass transfer efficiency (3. The initial primary mass is 5M for the figures 
on the left and 25M Q for the figures on the right. In the two top figures conservative evolution is assumed, /3 = 1.0, and the two 
lower figures f3 = 0.5. The subtypes abbreviations are described in the text. 



thermal equilibrium. This leads to a contact situation dur- 
ing the first rapid phase of case A mass transfer. If the mass 
transfer rate M just before contact occurs is larger than 2% 
of the typical thermal mass transfer rate M/tkh a system is 
classified as case AR. 

- Case AS Contact during slow mass transfer. Binaries 
with small periods and moderately equal masses stay semi- 
detached during thermal-timescale mass transfer, but come 
into contact during the subsequent slow nuclear-timescale 
mass transfer phase. 

- Case AO Contact by overtaking: Mass transfer accelerates 
the evolution of the mass-accreting star. Especially in lower- 
mass systems there is a significant range of initial separations 
for which the mass-accreting star depletes its central hydro- 
gen first. Most systems come into contact shortly afterwards 
as the mass-accreting star expands in t he Hertzsprung gap. 
This s ubtype corresponds to case AE in lNelson & Eggletonl 
(1200 lh . with the small difference that we explicitly check if 
the system comes into contact. 

- Case AN No contact on the MS: Wider systems with mod- 
erately equal mass can avoid contact while both components 
are still on the main sequence. This subtype comprises cases 



AG, AL, AN and AB in lNelson & Eggletonl (1200 lh between 
which we do not distinguish here. 



3.2. Comparison to \Nelson & Egqletori \2001\) 

In contrast to iNelson & Eggletonl feOOlb our models are calcu- 
lated at Z = 0.004 instead of Z = 0.02. Due to lower opacity the 
stars are more compact, which causes the different subtypes of 
case A binary evolution to shift to smaller orbital periods. 

Another diff er ence between our survey and 
INelson & Eggletonl d200ll) is that in our version of the 
code the equations for both stars are solved simultaneously 
instead of first evolving the primary star, keeping track of the 
mass transfer rate and then evolving the secondary applying the 
stored mass transfer rate. This enables us to model reverse mass 
transfer, although this does not occur often and it is not very 
relevant for the comparison with observations we undertake. A 
second advantage is that we continue our evolutionary calcula- 
tions when the stars reach contact, see Figure [2] Therefore we 
find a subset of type AR systems that come into contact during 
rapid mass transfer, but in which the mass-accreting star quickly 
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Fig. 2. The natural logarithm of the radius over the Roche lobe 
radius as function of time for the mass-accreting star in binary 
systems with initial primary mass Mij = 18M , initial orbital 
period f, = 2.2 days for three different mass ratios. The primary 
star of the system with logq = 0.2 fills its Roche lobe at an 
age of « 8.87 Myrs. During the phase of rapid mass transfer the 
mass-accreting star expands so much that contact is established. 
Shortly afterwards it restores its thermal equilibrium and shrinks 
inside its Roche lobe again. The mass-accreting star in the sys- 
tem with logg, =0.1 stays well inside its Roche-lobe. The ac- 
creting star in the log q, = 0.3 system excess its Roche Lobe 
by more than 10%. We stop the evolutionary calculations as our 
code is not able to handle deep contact configurations. Systems 
with more extreme initial mass ratios are likely to merge during 
rapid mass transfer. 

restores its thermal equilibrium and shrinks inside its Roche 
lobe again. This is followed by a slow mass-transfer phase 
like in systems that avoid contact altogether during rapid mass 
transfer. Some of the observed systems may have followed this 
type of evolution. Further evolution is either as case AS, AO or 
AN. 

3.3. The influence of the mass transfer efficiency 

An important effect of lowering the mass transfer efficiency 
is that more systems avoid contact during rapid mass transfer. 
The critical mass ratio, which defines the border between type 
AN and type AR systems is up to 25% larger if we assume 
B = 0.5 than for the conservative tracks, see Fig Q] Since con- 
tact during rapid mass transfer potentially has catastrophic con- 
sequences for the further evolution of the binary, we expect that 
non-conservative mass transfer results in a larger number of bi- 
naries reaching a long-lived semi-detached phase. 

The mass transfer efficiency also strongly influences how 
much time a binary systems spends in a semi-detached phase, 
see Fig [3] The accreting star will be less massive in non- 
conservative models and therefore it will evolve more slowly. 
The duration of the slow mass transfer phase can be twice as long 
for non-conservative models (B = 0.5) compared to the conser- 
vative models. 

The implication of non-conservative mass transfer that more 
binary systems become semi-detached and the duration of the 
semi-detached phase is longer, can in principle give a strong 
test of the mass transfer efficiency if the number of semi- 
detached systems is compared to the number of detached sys- 



tems. However, a good understanding of the selection effects 
and knowledge of the initial mass-ratio distribution are critical 
in performing such a test. 

4. Observed systems 

The stellar masses, radii and temperatures of the components 
of 50 double- lined eclip sing O B-type binaries in the SMC were 
presented by lH03l and IH05L In this section we will summa- 
rize their work. The systems were initially taken from 1400 
eclipsing binaries detected in the Optical Gravitational Lensing 
Experiment (OGLE) survey. 169 systems were chosen for spec- 
troscopic follow-up observations by selecting systems brighter 
than B < 16, to allow a sufficient signal-to-noise ratio, and 
orbital period P or b < 5, to ensure adequate phase coverage. 
About 100 systems had spectroscopic observations obtained 
near quadratu re phases calculated from the adopted photom et- 
ric ephemeris dWvrzvkowski et alJl2004tlUdalski et al.|[l998l) . A 
spectral disentangling procedure was used to establish the or- 
bital parameters directly, together with the average disentangled 
spectrum of each binary component. Orbital solutions were de- 
termined for a total of 50 systems, and were combined with 
analysis of the OGLE /-band light curves to yield complete as- 
trophysical parameters. The 50 systems are spread over the 10 
OGLE fields which cover the central 2.4 square degree area of 
the Small Magellanic Cloud. 

Twenty-one systems are in a detached configuration, two are 
contact systems. The remaining 27 systems are semi-detached 
and believed to be undergoing their first phase of mass trans- 
fer. In Table Q] the observed stellar parameters are summarized. 
For two of the detached systems, 1 and 6, alterna tive s emide- 
tached solutions were found after the publication of lH05l These 
are listed as system 51 and 52 respectively. 

The light curves show no evidence for departures from cir- 
cular orbits, except for the following detached systems: 10 (e = 
0.063+0.002), 14 (e = 0.113+0.004), 16 (e = 0.190+0.006) and 
19(e = 0.013+0.009). We can therefore assume that the spin pe- 
riods of the stars in each system are synchron ized with the orbital 
period . Th e orbital periods are taken from IWvrzvkowski et al.l 
d2004 and lUdalski etail dl998l) . The mass of the brightest com- 
ponent, the stellar radii, their errors and the temp erat ure of the 
brightest component are taken directly from |H03| and lH05l We 
derived the mass ratio q = M p /M s (where M p is the mass of the 
primary, i.e. the star eclipsed during the deepest eclipse which 
is usually the most luminous component, and M s is the mass 
of the companion) and its error directly from the radial velocity 
semi-amplitudes and K p and their errors, given in the same 
papers (q = K s /K p ). The three authors independently classified 
the spectral type of the primary component to determine the ef- 
fective temperature. They estimate the error on the temperature 
to be in the order of the difference in temperatur e betw een two 
spectral subtype^ An erro r of 1 000 K is given in lH03l while an 
error of 1500 K is given in lH05l However, the average temper- 
ature difference between spectral subtypes in the range 06-B2 
is cr Tclr ss 1800K. In addition there are significant uncertainties 
in the conversion of spectral type to temperature at the metal- 
licity of the SMC. Therefore we decided to adopt an error of 
Ar e ff = 2000 K to allow for some systematic un certai nty. The 
temperature ratio and its error were determined bvlHOlandlHOl 
from the /-band flux ratio. Unfortunately this band is not very 

3 The errors in the secondary temperatures quoted in lH03l and lH05l 
only reflect the uncertainty in the flux ratio and not the uncertainty in- 
duced by the error on the temperature of the primary star. 
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p = 1.0 M 1i = 5M p = 1.0 M 1i = 25M :l , 




0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.05 0.1 0.15 0.2 0.25 0.3 0.35 

Log (M v IM 2:i ) Log {M 1 JM z i ) 

Fig. 3. Time in years (color shading) spent in a semi-detached configuration for binary systems with different initial parameters. 
For a description of the difference between the four diagrams and their axis see Fig.Q] The duration of the semi-detached phase for 
non-conservative models is longer than for conservative models and more systems avoid contact during rapid mass transfer, see also 
Fig.Q] This trend is seen for all different primary masses in our grid of binary models over the full range of /3 



sensitive to temperature for OB stars as it covers the Raleigh- 
Jeans tail of the spectrum. Nevertheless, the formal errors on 
the temperature ratio as quoted in Table Q] are often very small. 
Photometric observations in additional bands or higher resolu- 
tion spectra would improve the determination of the primary 
temperature, the temperature ratio and their errors. 

5. Fit method 

For each system in the observed sample we determine the 
best fitting binary evolution track using a least-squares fit. Our 
four fitting parameters are the three initial binary parameters 
that determine the binary evolution track - the initial primary 
mass Mi j, the initial mass ratio q, and the initial orbital period P, 
- and the age t of the best fitting model on that track. 

There is freedom in the set of observed parameters we can 
use for fitting. To avoid propagation of errors, we preferentially 
use observables that are determined directly from the light curve 
and the spectra. For a proper ^- 2 -test the observables should be 
independent. In case of the semi-detached systems the radius of 
the Roche-lobe filling star is determined from the light curve, 
but it is not independent of the mass ratio and the orbital pe- 
riod. Therefore we cannot use all three observables at the same 
time. In detached binaries the orbital period does not change sig- 
nificantly during the evolution before the onset of mass trans- 
fer. Fitting the period within the tiny observed error is pointless 



given the much larger uncertainties in the other observables. We 
therefore choose a definition ofx 1 which does not include the or- 
bital period. We do demand that a fit solution matches the period 
within the accuracy of our grid. 

We use the following set of six observed parameters, listed 
in Table Q] the mass of the most luminous star logM p , the mass 
ratio log q, the radius of the most luminous star log R p , the radius 
of the companion star log7? s , the temperature of the most lumi- 
nous star log T p and the temperature ratio log r p /r s . With six in- 
dependent observables and four model parameters we have two 
degrees of freedom (v = 2) when we fit semi-detached systems. 
In the case of detached systems one of our model parameters, 
the initial orbital period, is degenerate because stars evolve es- 
sentially as if they were single until the onset of mass transfer. In 
practice we therefore have three model parameters which yields 
v = 3 for the detached systems. 

The cube of models only has discrete values for the model 
parameters M\ j, q, and P,. Ideally we want to vary the model pa- 
rameters continuously, but our grid spacing is not fine enough for 
some of the more accurate ly determined observables. Following 
iNelson & Eggletonl (l2001l) we therefore introduce 'theoretical 
errors' cr^ to account for the discreteness of the model grid. For 
the total error cry on the y'-th observed parameter we add the ob- 
servational error cr t, s and the model error <r± in quadrature, 
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Table 1. Parameters of 50 double-lined eclipsing binaries in the SMC. Details are given in Section|4] 
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id OGLE-id | M P (M Q ) M p /M s 


Rp(Re) 


RARe) 


T„(kK) 


T P IT S 


P(d) 


Detached systems 



01 


04 


056804 


13.0 


±0.3 


0.912 ±0.019 


4.11 


±0.09 


4.94 


±0.08 


30.1 ±2.0 


1.245 ±0.017 


1.08987 


02 


04 


103706 


17.5 


±0.6 


1.768 ± 0.020 


5.36 


±0.13 


4.05 


±0.12 


27.8 ± 2.0 


1.028 ±0.017 


1.35585 


03 


04 


163552 


13.3 


± 1.0 


1.070 ±0.052 


5.31 


±0.15 


5.18 


±0.15 


25.5 ± 2.0 


1.005 ±0.007 


1.54581 


04 


05 


038089 


17.1 


± 1.5 


0.893 ± 0.058 


6.27 


±0.21 


6.19 


±0.26 


30.1 ±2.0 


1.033 ±0.036 


2.38946 


05 


05 


095194 


20.3 


±4.5 


0.870 ±0.149 


8.06 


±0.56 


9.53 


±0.66 


33.8 ±2.0 


1.038 ±0.002 


3.18742 


06 


05 


140701 


6.9 


±0.7 


1.307 ± 0.064 


7.88 


±0.31 


7.01 


±0.31 


23.5 ± 2.0 


1.535 ±0.025 


3.62544 


07 


05 


180064 


10.7 


±0.4 


1.521 ±0.019 


5.57 


±0.27 


4.49 


±0.29 


25.5 ± 2.0 


1.493 ±0.107 


2.51491 


08 


05 


255984 


11.6 


±2.0 


1.659 ±0.133 


4.16 


±0.43 


3.38 


±0.34 


25.5 ± 2.0 


1.019 ±0.038 


1.56384 


09 


05 


305884 


17.6 


± 1.3 


1.086 ±0.059 


7.76 


±0.27 


6.52 


±0.37 


33.8 ±2.0 


1.040 ±0.007 


2.17648 


10 


05 


311566 


12.9 


±0.6 


1.261 ± 0.034 


4.38 


±0.13 


3.21 


±0.19 


30.1 ±2.0 


1.050 ±0.024 


3.29139 


11 


06 


011141 


15.1 


±0.3 


1.068 ± 0.016 


5.05 


±0.05 


5.01 


±0.06 


30.1 ±2.0 


1.324 ±0.027 


1.17737 


12 


06 


180084 


16.2 


± 1.3 


1.169 ±0.073 


5.85 


±0.25 


5.70 


±0.25 


27.8 ± 2.0 


0.995 ±0.011 


1.47523 


13 


06 


215965 


16.0 


±0.5 


0.931 ±0.027 


9.91 


±0.19 


10.35 


±0.19 


27.8 ± 2.0 


1.000 ±0.037 


3.94304 


14 


06 


221543 


11.9 


± 1.8 


1.025 ±0.103 


5.27 


±0.48 


4.53 


±0.52 


25.5 ± 2.0 


0.959 ± 0.009 


3.41678 


15 


07 


120044 


12.5 


±0.5 


1.004 ± 0.028 


4.78 


±0.10 


4.68 


±0.12 


25.5 ± 2.0 


0.998 ± 0.016 


1.31081 


16 


07 


255621 


9.3 


±0.8 


1.293 ± 0.049 


5.04 


±0.29 


3.54 


±0.56 


25.5 ±2.0 


1.091 ±0.043 


4.33063 


17 


08 


087175 


12.0 


± 1.0 


1.119 ±0.073 


4.48 


±0.21 


4.38 


±0.22 


25.5 ± 2.0 


1.033 ±0.021 


1.10226 


18 


08 


104222 


13.1 


±0.9 


1.078 ± 0.055 


5.25 


±0.38 


5.19 


±0.38 


25.5 ± 2.0 


1.012 ±0.033 


1.56384 


19 


10 


037156 


19.5 


±0.4 


1.147 ±0.014 


7.19 


±0.16 


4.41 


±0.46 


32.2 ± 2.0 


1.028 ±0.012 


2.69834 


20 


10 


110440 


10.8 


±2.1 


1.825 ±0.131 


3.97 ± 0.39 


3.66 


±0.34 


25.5 ± 2.0 


1.194 ±0.027 


1.56384 


21 


11 


057855 


12.4 


± 1.1 


1.510 ±0.061 


5.22 


±0.21 


3.73 


±0.24 


27.8 ± 2.0 


1.072 ±0.039 


1.29695 



Semi-detached systems and contact systems (47 and 49) 



22 


01 099121 


11.3 


±0.8 


1.704 ± 0.048 


4.97 


±0.19 


6.66 


±0.24 


27.8 ± 2.0 


1.758 ±0.018 


2.45890 


23 


04 110409 


13.7 


±0.8 


1.544 ± 0.046 


4.34 


±0.20 


8.36 


±0.29 


25.5 ± 2.0 


1.660 ±0.017 


2.97315 


24 


05 026631 


11.5 


±0.6 


1.021 ± 0.039 


5.12 


±0.11 


5.65 


±0.12 


25.5 ± 2.0 


1.489 ±0.026 


1.41169 


25 


05 060548 


10.8 


±0.4 


1.238 ± 0.027 


8.39 


±0.14 


9.63 


±0.15 


30.1 ±2.0 


1.722 ±0.025 


3.63863 


26 


05 202153 


19.9 


± 1.1 


1.592 ± 0.036 


9.53 


±0.26 


12.84 


±0.32 


32.2 ± 2.0 


1.374 ±0.012 


4.60677 


27 


05 208049 


10.0 


±0.2 


2.093 ± 0.014 


6.62 


±0.08 


6.82 


±0.08 


25.5 ± 2.0 


1.791 ±0.020 


3.02982 


28 


05 243188 


27.3 


± 1.5 


1.466 ± 0.041 


7.29 


±0.19 


7.89 


±0.20 


35.5 ±2.0 


1.104 ±0.002 


1.87174 


29 


05 277080 


17.4 


±0.9 


1.540 ±0.036 


5.05 


±0.14 


6.82 


±0.17 


25.5 ± 2.0 


1.607 ±0.015 


1.93934 


30 


05 300549 


25.4 


± 1.0 


1.458 ± 0.027 


6.35 


±0.11 


6.15 


±0.11 


30.1 ±2.0 


1.879 ±0.030 


1.33167 


31 


05 316725 


16.9 


± 1.5 


1.916 ±0.048 


6.13 


±0.25 


7.66 


±0.29 


33.8 ±2.0 


1.072 ±0.006 


2.55606 


32 


06 077224 


15.9 


± 1.0 


1.218 ± 0.049 


8.93 


±0.26 


11.52 


±0.32 


32.2 ± 2.0 


1.253 ±0.026 


3.82087 


33 


06 152981 


12.5 


±0.4 


1.533 ±0.018 


4.94 


±0.08 


6.26 


±0.09 


25.5 ± 2.0 


1.268 ±0.009 


2.00334 


34 


06 158118 


16.0 


±0.7 


2.029 ± 0.023 


7.61 


±0.15 


7.38 


±0.15 


29.2 ± 2.0 


1.589 ±0.020 


2.57832 


35 


06 251047 


8.1 


±0.2 


1.473 ± 0.016 


4.62 


±0.11 


6.40 


±0.12 


23.5 ±2.0 


1.738 ±0.012 


2.51323 


36 


06 311225 


21.2 


±0.4 


1.773 ±0.012 


6.55 


±0.06 


6.69 


±0.06 


30.1 ±2.0 


1.282 ±0.015 


1.84606 


37 


06 319960 


10.6 


±0.8 


1.569 ± 0.057 


4.47 


±0.18 


9.38 


±0.34 


25.5 ± 2.0 


1.567 ±0.014 


4.05690 


38 


07 066175 


19.6 


± 1.8 


1.703 ± 0.059 


7.88 


±0.31 


10.38 


±0.40 


32.2 ± 2.0 


1.271 ±0.021 


3.62638 


39 


07 142073 


12.6 


± 1.2 


2.000 ± 0.053 


9.57 


±0.35 


7.75 


±0.28 


30.1 ±2.0 


1.435 ±0.019 


3.18871 


40 


07 189660 


15.3 


± 1.2 


1.498 ± 0.046 


5.35 


±0.16 


5.98 


±0.18 


25.5 ± 2.0 


1.510 ±0.014 


1.66784 


41 


07 193779 


11.6 


± 1.0 


1.966 ± 0.053 


5.79 


±0.21 


4.93 


±0.18 


25.5 ± 2.0 


1.486 ±0.019 


1.67293 


42 


07 243913 


18.6 


± 1.1 


1.768 ± 0.034 


8.23 


±0.20 


8.28 


±0.20 


32.2 ± 2.0 


1.242 ±0.015 


2.63160 


43 


08 209964 


18.8 


±0.9 


1.300 ±0.031 


7.67 


±0.20 


10.69 


±0.24 


36.3 ± 2.0 


1.371 ±0.014 


3.31475 


44 


09 010098 


17.8 


± 1.8 


1.304 ± 0.075 


5.08 


±0.22 


5.06 


±0.22 


33.8 ±2.0 


1.062 ±0.003 


1.11208 


45 


09 047454 


12.6 


±2.3 


1.362 ±0.107 


4.20 


±0.26 


5.58 


±0.35 


27.8 ± 2.0 


1.644 ±0.029 


1.57397 


46 


09 064498 


8.4 


±0.7 


3.092 ± 0.044 


5.40 


±0.24 


5.10 


±0.23 


25.5 ± 2.0 


1.496 ±0.032 


2.63512 


47 


09 175323 


23.5 


± 1.6 


1.451 ±0.041 


10.16 


±0.26 


8.48 


±0.23 


39.3 ± 2.0 


1.021 ±0.006 


2.20596 


48 


10 094559 


12.0 


± 1.0 


1.208 ± 0.065 


4.68 


±0.21 


6.19 


±0.26 


30.1 ±2.0 


1.262 ±0.019 


1.75193 


49 


10 108086 


16.9 


± 1.2 


1.183 ±0.057 


5.67 


±0.18 


5.32 


±0.16 


33.8 ±2.0 


1.130 ±0.038 


0.88308 


50 


11 030116 


14.3 


± 1.9 


1.868 ± 0.086 


7.42 


±0.42 


8.05 


±0.46 


27.8 ± 2.0 


1.368 ±0.016 


2.95427 



Alternative solutions for system 1 (51) and system 6 (52) 



51 


04 056804 


16.8 


±2.1 


1.724 ± 0.088 


4.46 


±0.25 


4.48 


±0.25 


30.1 


±2.0 


1.250 ±0.009 


1.08987 


52 


05 140701 


7.0 


±0.7 


1.307 ± 0.064 


6.00 


±0.20 


8.30 


±0.30 


24.3 


±2.0 


1.493 ±0.007 


3.62544 



This leads to the following definition of^f 2 , 

o Xobs ; - X th i ( Mi U q„ Pi, t) ] 2 

X 2 (M u , q u P„ t) = - V J ■ (9) 

v *-H cr 1 . 
j i 

The sum is taken over all observed parameters X \, s j and the 
corresponding values in our models X t hj. We determine the best 
fitting evolution track (Mj ., q*, PI) and the best fitting age (f) 



by minimizing x 2 with a systematic grid search of the fitting 
parameters. We define 

xL=x*(M* 1J ,qir i .n do) 

We only allow fit solution in which the observed period is 
matched within |logP t, s - logP t h| < 0.05. For the theoretical 
errors we take half the initial grid spacing, i.e. for logMi i and 
log Pj, cr th = 0.025 and for \ogq„ cr th = 0.0125. We estimate 
the initial grid spacing for the radii and the effective temperature 
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for each system individually from the two initial models with 
primary masses closest to the observed mass of the primary. 

The initial grid spacing for the logarithm of the temperature 
ratio scales approximately linearly with the spacing of the loga- 
rithm of the mass ratio, with no dependence on the orbital period 
and a weak dependence on the primary mass. For a typical pri- 
mary mass of 1OM we fit a straight line through log T1/T2 as 
function of log M\ /M2, which results in an estimate for the theo- 
retical error on the logarithm of the temperature ratio of 0.0065. 

We realize that the introduction of theoretical errors in the 
definition ofx 2 compromises a rigorous statistical interpretation 
of the obtained v 2 . values. Furthermore, using the initial grid 
spacing is not ideal, especially for post-mass transfer systems, 
because the real grid spacing will vary during the evolution. 
Therefore we compared different definitions of x 2 and of <r±. 
We tested the effect of using the orbital period as an observable 
instead of the radius of the Roche-lobe filling star, and we varied 
the magnitude of cr th between zero a nd the full initial grid spac - 
ing. The latter definition was used bv lNelson & Eggle tonl(l2001l) . 
The magnitude x\^ m ^ or a certam system depends directly the 
chosen definition; in particular, overestimating cr t h leads to arti- 
ficially low values of x 2 - On the other hand, o-± = can yield 
artificially high values of x 1 due to incomplete sampling of the 
model parameters. Therefore, although we attempted to choose a 
reasonable definition, x 2 . can not be used as an objective crite- 

' 1 mm J 

rion to determine whether a fit is good or bad, only as a relative 
estimate. Nevertheless, we found that the y 2 . values obtained 

' a mm 

with different methods are well correlated; in other words, the 
relative quality of a fit is independent of the chosen fit method. 
The best fitting parameters are not sensitive to the chosen defi- 
nition, i.e. do not vary more than one step in our grid of models 
except for a few systems for which no good fit could be obtained 
at all. 

For the detached systems we only take detached solutions, 
i.e. the radii of both stars should be smaller than 98% of their 
Roche-lobe radii. This definition i s con sistent with the defini- 
tion used for the observed systems (|H05|) . For the semi-detached 
systems we only accept solutions during which the evolution is 
slow, on a nuclear timescale. Case A mass transfer starts with 
a phase of rapid mass transfer on the thermal timescale, Tkh- 
After the mass ratio reverses, the orbit widens and a phase of 
slow mass transfer follows on the nuclear timescale, r nuc , of the 
mass-donating star. As Tkh ~ 0.01t 11uc one would expect only 1 
out of 100 systems to be in the rapid mass transfer phase, if no 
other biases play a role. We therefore assume that all observed 
semi-detached systems are currently in the slow mass transfer 
phase, and we only fit models in which: 

- the donor star is filling its Roche lobe, such that the radius 
of the initially more massive star is larger than 98% of the 
Roche-lobe, 

- the stars are in thermal equilibrium, i.e. the thermal luminos- 
ity of the mass donor is less than 1 % of its total luminosity, 

- the mass transfer rate is smaller than 1% of the typical 
thermal-timescale mass transfer rate of the donor, Md/t^kh, 
and 

- the donor star has become less massive than the accreting 
star, M D < M A . 

To determine confidence limits on the fitted model parame- 
ters one has to determine for which values of the parameters the 
condition x 2 (Mi, qi, P,, f) < x 2 m { n + c is satisfied, varying each 
parameter in turn, where c — 1/v to determine lcr boundaries 
for a propers-test. To compensate for any systematic over- or 



underestimation of the observational and theoretical errors we 
take c — x 1 ■ /v. As our model grid is discrete we determine the 

^ min ' fc 

upper and lower limit on each parameter by the maximum and 
minimum values for which the condition is satisfied, plus or mi- 
nus half of the grid spacing for that parameter. To estimate the 
range of f3 values which fit well for a certain model, we first de- 
termine the best fitting model for each of the four assumptions 
of /3, and the corresponding X^m- T° °^ tam ^min as a function of 
P on the full range 0.25 < /3 < 1.0 we interpolate linearly be- 
tween the four values. Let /3* denote the best-fitting value of /?. 
To estimate the confidence limits on /3 we determine for which 
values x 1 ■ (P) < X 2 ■ (fi*) + Y 2 - (J3*)/v. 

6. Results 

The best fitting model parameters for the detached systems are 
given in Table [2] and for the semi-detached systems in Table [3] 

6.1. Detached Systems 

For twelve of the twenty-one detached systems we find good or 
reasonable fit solutions with < 3, see Table [2] Inspection 

by eye confirms that the fits for these systems are indeed reason- 
ably good, see Figure|4]for an example. An additional set of four 
systems have 3 < x 2 ■ < 5, which we still regard as acceptable 
fits on the basis of eye inspection. For five systems no acceptable 
fit could be obtained C^ in > 10). In Sect. l6.1.Tl we comment on 
the individual fit solutions of these systems. Figures similar to 
Fig. |4]for all systems are given in the on-line Appendix. Here 
we first discuss the overall properties of the results as presented 
in Table [2] 

If we eliminate system 1 because of the existence of a well 
fitting alternative semi-detached solution and system 7 because 
of bad quality of the data (see Sect. 16.1. It we are left with good 
or acceptable fits for 16 out of 19 (84%) of the detached systems. 
We may thus conclude that on the whole our models are able to 
reproduce the detached binaries in the sample quite well. 

In some cases the best fitting solution lies on the border of 
our grid. This happens often for systems with a mass ratio near 
1 as our grid only holds models with initial mass ratios down 
to 1.059. In two cases the observed orbital period is too high 
for our grid. This is the case for system 6, for which no good fit 
could be obtained for other reasons than the period alone, and for 
system 16 which nevertheless gives a good fit. The boundaries 
for the masses, the lower boundary for the period and the upper 
boundary for the mass ratio are never reached. 

Out of the 16 systems for which we obtained good or accept- 
able fits three systems are so wide that our conservative models 
predict that mass transfer will not start before the most massive 
star enters the Hertzsprung gap, six systems will start mass trans- 
fer on the main sequence but avoid contact, at least until the pri- 
mary reaches the Hertzsprung gap, and six systems will come 
into contact during slow mass transfer, see the last column of 
Table [2] One system will reach contact already during the rapid 
mass transfer phase. 

6.1 .1 . Comments on individual fits 

We first give brief notes on the systems that fit well and then 
more detailed notes on the systems for which no good fits could 
be obtained. 

Several systems have lower observed temperatures than our 
models predict: 3, 4, 15 and 18 and to a lesser extent 12 and 17. 
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Fig. 4. Example of a fit for a detached system. The star-shaped symbols indicate the observed parameters for system 21, the 
circles (green in on-line version) indicate the best-fitting model. The conservative evolution track of the best-fitting model is plotted 
in dark gray (red) for the star indicated as the primary, see Table [T] and in lighter gray (blue) for the companion. On the left a 
Hertzsprung-Russel diagram is given for reference, although we do not use the luminosity as a fit parameter. The error bars on the 
filled star symbols indicate the error on the primary mass, temperature and radius. In contrast the error bars on the open star symbols 
indicate the error on the mass ratio, the temperature ratio and the secondary radius. The error bars on the circles indicate the adopted 
theoretical errors on the same quantities (see Sect. I5J, which correspond to half the initial grid spacing of the models. 



Table 2. Fit-results for the detached systems: The best fitting initial primary mass My, initial mass ratio q\ - Mx^jMi^ and initial 
orbital period P\ are given together with estimates of the lcr confidence limits (for details see the last paragraph of Section|5]l. The 
last two columns give the best fitting age on the best fitting evolution track and, for systems with^f 2 < 10, the evolutionary subtype 
assuming conservative mass transfer, see Section [3~T1 The notation AR — » ... indicates that the best-fitting model experiences a 
temporary contact phase during rapid mass transfer after which the evolution continues as indicated. 



id 


2 


M u (M ) 




P, (days) 


/ (Myrs) 


case 


1 


43.0 


14.1 


-4.7 


+4.7 


1.05 


-0.05 


+0.24 


1.12 


-0.21 


+0.18 


3.1 




2 


21.8 


15.8 


-2.5 


+5.3 


1.58 


-0.21 


+0.15 


1.31 


-0.16 


+0.32 


6.8 




3 


1.2 


12.6 


-0.7 


+0.7 


1.06 


-0.06 


+0.03 


1.51 


-0.08 


+0.29 


10.6 


AS 


4 


1.7 


17.8 


-2.8 


+ 1.0 


1.06 


-0.06 


+0.03 


2.18 


-0.12 


+0.60 


6.2 


AN 


5 


1.0 


20.0 


-1.2 


+3.7 


1.13 


-0.13 


+0.02 


3.21 


-0.39 


+0.60 


8.1 


AN 


6 


114.0 


10.0 


-3.3 


+27.6 


1.06 


-0.06 


+ 1.31 


3.47 


-0.42 


+ 1.20 


23.7 




7 


10.7 


11.2 


-1.8 


+3.8 


1.49 


-0.12 


+0.14 


2.57 


-0.41 


+0.59 


14.8 




8 


4.5 


10.0 


-1.6 


+3.3 


1.33 


-0.11 


+0.21 


1.55 


-0.18 


+0.40 


15.3 


AN 


9 


1.2 


20.0 


-3.2 


+ 1.1 


1.13 


-0.04 


+0.02 


2.27 


-0.36 


+0.47 


6.3 


AN 


10 


2.3 


12.6 


-0.7 


+0.7 


1.26 


-0.11 


+0.04 


3.02 


-0.17 


+0.78 


3.9 


AN 


11 


30.9 


15.8 


-3.9 


+5.3 


1.18 


-0.18 


+0.19 


1.31 


-0.30 


+0.20 


4.3 




12 


2.0 


14.1 


-0.8 


+2.7 


1.05 


-0.05 


+0.10 


1.58 


-0.25 


+0.27 


9.4 


AS 


13 


3.6 


17.8 


-2.8 


+ 1.0 


1.06 


-0.06 


+0.03 


4.34 


-0.84 


+0.53 


10.0 


B 


14 


2.1 


11.2 


-1.8 


+2.1 


1.06 


-0.06 


+0.10 


3.63 


-0.71 


+0.63 


12.7 


B 


15 


2.0 


12.6 


-2.0 


+0.7 


1.06 


-0.06 


+0.03 


1.20 


-0.07 


+0.42 


8.0 


AS 


16 


0.3 


10.0 


-0.6 


+0.6 


1.27 


-0.04 


+0.11 


3.89 


-0.22 


+0.48 


17.1 


B 


17 


0.9 


11.2 


-0.6 


+0.7 


1.06 


-0.06 


+0.10 


1.15 


-0.18 


+0.14 


9.7 


AS 


18 


0.9 


12.6 


-0.7 


+0.7 


1.06 


-0.06 


+0.03 


1.51 


-0.08 


+0.39 


10.7 


AS 


19 


3.1 


17.8 


-1.0 


+3.3 


1.19 


-0.10 


+0.03 


2.44 


-0.14 


+0.77 


7.2 


AN 


20 


3.7 


11.2 


-1.8 


+2.1 


1.58 


-0.12 


+0.15 


1.45 


-0.08 


+0.51 


9.6 


AR 


21 


1.4 


12.6 


-0.7 


+0.7 


1.42 


-0.04 


+0.12 


1.35 


-0.21 


+0.16 


10.0 


AR^ AS 



Good fits with small values of v 2 . were nevertheless obtained, 

A min ' 

because of the large adopted error on the effective temperature 
of the primary. However, the discrepancy is systematic and the 
opposite situation does not occur. In addition, the components 
of system 12 have very equal temperatures, while the masses are 
different. 

Both s tars o f system 10 are very close to the zero-age main 
sequence. HQ5 note that the orbit of this system is slightly eccen- 
tric, see Section [4] which is consistent with the young age and 
the fact that both stars are well inside their Roche lobes. Also 
systems 14 and 16 have significant eccentricities. Although the 



stars in these binaries are more evolved than in system 10, both 
systems have wide orbits and are well detached. 

For four systems the fits we obtained were only marginally 
acceptable, 3 <;t min — 5: 

8 For this system we cannot fit the small temperature ratio well 
at the same time as the moderately large mass ratio, resulting 
iny 2 . =4.5. 

A mm 

13 In this wide and evolved binary, the equal temperatures and 
the large and almost equal radii of the stars suggest that the 
masses must also be nearly equal. The best fit (x^m ~ ^.6) 
is obtained with the smallest mass ratio available in our grid, 
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q — 1.06. We expect that if a model with equal masses were 
available, we would have found a very good fit for this sys- 
tem. 

19 The stars in this system have fairly equal masses, while the 
radii are very different. The less mass ive co mponent has the 
size of a zero-age main-sequence star. |H05| note that the sys- 
tem may be slightly eccentric. One possibility is that the 
larger star is still a pre-main sequence star, but this is also 
the more massive star which should have contracted towards 
the zero-age main sequence first. Another speculative possi- 
bility is that the stars were not born at the same time and 
that they became part of a binary system after the larger 
star started burning hydrogen. However, given the moderate 
value ofx 2 - - 3.1 such exotic scenarios are probably not 
required. 

20 This system has a large secondary radius compared to its 
mass (the radii are nearly equal whereas the masses are not) 
which cannot be explained by our detached models. As a re- 
sult the overall fit is not very good Otf^ in = 3.7). An alterna- 
tive fit solution was found for this system during a post-mass 
transfer phase. After a rapid phase of mass transfer the donor 
star briefly detaches when it restores its thermal equilibrium. 
This alternative fit was not significantly better and the prob- 
ability of observing a system in the short post-mass transfer 
phase is small compared to the pre-mass transfer solution. 

For five systems no good fit could be obtained: 

1 This system has nearly equal masses, which is not consistent 
with the large temperature ratio and the very different radii. 
The brightest star is the less massive component of the sys- 
tem which fits well to the evolution track of a star that has 
just arrived on the main sequence. However, the less bright 
and more massive component is too big and too cool com- 
pared to our models. 

For this system an alternative semi-detached light curve so- 
lution was found. This alternative solution can be fitted well 
against the semi-detached models with a reduced;^ 2 ^ o f 1 .6. 
We propo se to reject the detached solution published in lH03l 
and lH05l and adopt the new semi-detached solution indicated 
as system 51, see Sect. 16.21 

2 The mass ratio of the components in this system is large 
(M p /M s = 1.8) while the /-band flux ratio indicates that the 
effective temperatures of both stars are very similar. The less 
massive component is consistent with a single star model, 
while the more massive component is too cool and too small 
for its mass compared to the models. 

6 The large radii of both components in this systems suggest 
that the stars are close to hydrogen exhaustion in the core, at 
the end of the main sequence or even at the beginning of the 
Hertzsprung gap. This implies that, if no mass transfer has 
taken place yet, the masses of the stars must be nearly equal, 
approximately 6M . However, one would expect approxi- 
mately equal effective temperatures for both stars, which is 
not the case. An alternative semi-detached solution is avail- 
able for this system, but no acceptable fit against semi- 
detached models was obtained either (see Section |6\2l . 

7 The more massive component of this system can be fitted 
to the evolution track of a 1 1 M star, but the less massive 
component is cooler and bigger than one would expect for 
its mass. However, the light curve has a very small amplitude 
(0.06 mag) and is poorly defined due to large scatter (H05). 
Therefore the light curve solution is not very reliable and 
this system probably should not have been included in the 
sample. 



1 1 The equal masses of this system together with the extreme 
temperature ratio make this system difficult to fit. The more 
massive component is consistent with a single star of about 
15 M , but the less massive star is too cool compared to the 
models. 



6.2. Semi-detached Systems 

In total there are 29 semi-detached systems and 2 contact sys- 
tems in the sample, including the two detached systems (system 
1 and 6) for which we have alternative semi-detached solutions. 
The best fitting parameters for these systems are listed in Table[3] 
An example of a semi-detached system for which we obtained an 
excellent fit is given in Figure[5] 

Our models provide a good or reasonable fit to 11 out of the 
31 systems, with x\& n < 3, for at least one value of This in- 
cludes the two contact binaries. Four systems have 3 < x 2 . < 5, 
and another two have v 2 . < 6.5, and we consider these fits as 

^ nun 

(marginally) acceptable. The remaining 14 semi-detached sys- 
tems cannot be fitted well (y 2 . > 10). As we have mentioned 

^min ' 

before the absolute value of ^ /2 nin is not a very good measure of 
the goodness-of-fit. The distinction we make here, and the de- 
cision to accept fits with^ 2 iin < 6.5, is primarily guided by eye 
inspection of the solutions. In Sections |6.2.1| and |6.2.2| we briefly 
discuss the individual fit solutions, for good or acceptable and 
bad fits, respectively. 

There are several common problems that occur in fitting our 
models to semi-detached binaries. We distinguish three types of 
conflicts between the models and the observations: 

- q-conflict: The most ubiquitous discrepancy, occurring in 
more or less extreme form in all systems with x^ in > 5 ex- 
cept one, is the combination of a relatively small mass ratio 
(close to unity) with a large or very large temperature ratio. 
Such a combination cannot be achieved by our models for 
any value of /3. We discuss possible causes for this problem 
in Sect. 

- R-conflict: A second discrepancy that plagues nearly half 
the systems with^ 2 iin > 10, is that the ratio of donor radius 
over accretor radius, Ru/Ra, is smaller than shown by our 
models, in particular the radius of the accreting star is often 
under-predicted. 

- MT-conflict: A third problem affecting several systems, is 
the combination of a large primary mass with a relatively 
low temperature which cannot be matched simultaneously 
by our models. The inverse (small mass together with high 
temperature) also occurs in a few systems. 

We refer to these conflicts when discussing the individual sys- 
tems in Sect. 16.2.21 

The effect of varying the mass transfer efficiency on the fit 
solutions is often rather obscure, but we can make a number of 
general statements. The most obvious effect is that lowering ft 
results in a smaller (less extreme) post-mass transfer mass ra- 
tio, for given initial masses and period and at a given age. The 
same is true, although somewhat less obvious, for the temper- 
ature ratio. Since we are trying to fit a certain observed mass 
(and temperature) ratio, this means that models with smaller [5 
require either a larger (more extreme) initial mass ratio, or a 
larger amount of mass to be transferred, i.e. a more advanced 
evolution stage. A consequence of the latter is in many cases a 
larger donor radius compared to the radius of the accretor. This 
effect is often (though not always) seen in our best-fitting solu- 
tions for a certain system: a smaller j3 leads to a larger ratio of 
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Fig. 5. Example of a fit for a semi-detached system. The observed parameters for system 50 are plotted together with the best fitting 
model assuming a mass transfer efficiency of 0.5. See Fig|4]for description of the symbols. According to the best fitting model this 
system consisted initially of two stars with masses of about 13 and 12 M and an orbital period of 1.7 days. The most massive star, 
which lost about 5 M by Roche-lobe overflow, is now close to the end of its main sequence. This system will come into contact 
when the mass donor enters the Hertzsprung gap. 

Table 3. Fit-results for semi-detached and contact systems. See the caption of Table[2]for a description. 



id 


x l - 


P 


M u (M e ) 


If 


Pi (days) 


t (Myrs) 


case 


22 


16.9 


1.00 


7.9 


-1.3 


+1.5 


1 


06 


-0.06 


+0.16 


1.99 


-0.32 


+0.22 


33.9 




23 


17.4 


0.50 


10.0 


-1.6 


+ 1.9 


1 


19 


-0.19 


+0.26 


2.46 


-0.48 


+0.27 


23.2 




24 


50.0 


0.75 


10.0 


-3.3 


+5.0 


1 


12 


-0.12 


+0.33 


1.74 


-0.58 


+0.43 


20.6 




25 


97.5 


0.25 


11.2 


-2.8 


+7.6 


1 


06 


-0.06 


+0.24 


2.89 


-1.09 


+ 1.04 


19.7 




26 


4.1 


0.50 


17.8 


-1.0 


+3.3 


1 


06 


-0.06 


+0.03 


3.45 


-0.42 


+0.36 


10.2 


AN 


27 


26.4 


0.50 


7.1 


-0.4 


+2.4 


1 


05 


-0.05 


+0.10 


1.34 


-0.16 


+0.87 


50.1 




28 


1.8 


0.75 


25.1 


-1.4 


+4.8 


1 


33 


-0.17 


+0.13 


1.75 


-0.28 


+0.18 


5.6 


AN 


29 


31.1 


1.00 


10.0 


-1.6 


+3.3 


1 


12 


-0.12 


+0.25 


1.95 


-0.39 


+0.59 


21.6 




30 


141.5 


1.00 


11.2 


-3.7 


+7.6 


1 


06 


-0.06 


+0.57 


2.04 


-0.96 


+1.08 


18.2 




31 


1.8 


0.50 


20.0 


-1.2 


+ 1.1 


1 


69 


-0.06 


+0.03 


2.02 


-0.11 


+0.12 


10.7 


AR^ AN 


32 


5.7 


0.25 


17.8 


-1.0 


+3.3 


1 


06 


-0.06 


+0.10 


3.07 


-0.49 


+0.32 


9.9 


AN 


33 


1.4 


0.50 


12.6 


-0.7 


+0.7 


1 


26 


-0.04 


+0.04 


1.51 


-0.08 


+0.09 


14.7 


AN 


34 


15.3 


1.00 


12.6 


-2.0 


+2.4 


1 


06 


-0.06 


+0.10 


1.91 


-0.48 


+0.36 


15.8 




35 


41.4 


0.50 


7.1 


-1.8 


+2.4 


1 


12 


-0.12 


+0.18 


1.90 


-0.38 


+0.47 


42.7 




36 


3.7 


1.00 


15.8 


-0.8 


+3.0 


1 


11 


-0.11 


+0.26 


1.48 


-0.08 


+0.28 


9.9 


AS 


37 


3.1 


0.25 


10.0 


-0.6 


+0.6 


1 


19 


-0.04 


+0.11 


2.46 


-0.14 


+0.15 


23.7 


AN 


38 


0.7 


0.25 


17.8 


-1.0 


+ 1.0 


1 


13 


-0.04 


+0.03 


2.18 


-0.12 


+0.13 


10.5 


AN 


39 


10.4 


0.25 


11.2 


-0.6 


+2.1 


1 


06 


-0.06 


+0.03 


1.45 


-0.08 


+0.27 


22.7 




40 


17.7 


1.00 


10.0 


-2.5 


+ 1.9 


1 


06 


-0.06 


+0.23 


1.95 


-0.56 


+0.21 


21.6 




41 


3.0 


1.00 


8.9 


-1.4 


+0.5 


1 


06 


-0.06 


+0.09 


1.32 


-0.21 


+0.08 


26.9 


AO 


42 


0.6 


0.50 


17.8 


-1.0 


+ 1.0 


1 


06 


-0.06 


+0.03 


1.73 


-0.10 


+0.10 


9.6 


AN 


43 


14.1 


0.50 


17.8 


-2.8 


+3.3 


1 


13 


-0.13 


+0.17 


2.74 


-0.26 


+0.66 


9.7 




44 


1.2 


1.00 


17.8 


-2.8 


+ 1.0 


1 


19 


-0.19 


+0.18 


1.09 


-0.17 


+0.06 


3.4 


AS 


45 


11.6 


1.00 


7.9 


-2.0 


+2.7 


1 


12 


-0.12 


+0.33 


1.77 


-0.44 


+0.29 


32.7 




46 


6.5 


1.00 


7.1 


-0.4 


+ 1.3 


1 


68 


-0.22 


+0.05 


1.34 


-0.21 


+0.15 


63.5 


AR^ AN 


47 


2.4 


0.25 


28.2 


-1.6 


+5.3 


1 


19 


-0.10 


+0.11 


1.44 


-0.08 


+0.15 


6.1 


AS 


48 


3.5 


0.50 


11.2 


-1.8 


+2.1 


1 


18 


-0.09 


+0.19 


1.82 


-0.39 


+0.20 


17.2 


AN 


49 


0.9 


1.00 


15.8 


-0.8 


+ 1.0 


1 


05 


-0.05 


+0.04 


1.17 


-0.07 


+0.07 


5.9 


AS 


50 


1.2 


0.50 


12.6 


-0.7 


+2.4 


1 


06 


-0.06 


+0.03 


1.70 


-0.09 


+0.18 


16.3 


AN 


51 


2.6 


1.00 


12.6 


-2.0 


+0.7 


1 


42 


-0.19 


+0.04 


1.07 


-0.11 


+0.36 


12.0 


AR^ AS 


52 


14.0 


0.25 


8.9 


-1.4 


+ 1.7 


1 


06 


-0.06 


+0.09 


2.34 


-0.46 


+0.26 


28.4 





donor to accretor radius, as well as (or sometimes instead of) a 
less extreme mass ratio. 

A few systems can be fitted well for different assumptions of 
j6, such that we cannot determine, on basis of the available ob- 
servables, if the system has gone through conservative or non- 
conservative mass transfer. However, in some of these cases (i.e. 
systems 31, 33, 38 and 51) the surface abundances of carbon and 
nitrogen differ substantially between the best-fitting models for 



different [5. This is a consequence of the fact that a smaller fl 
yields a fit at a more advanced evolution stage, when the donor 
star is peeled off to a larger extent. Abundance determinations 
of carbon and nitrogen would therefore enable us in principle 
to constrain the mass transfer efficiency. Unfortunately these are 
not available at this moment, because the resolution of the spec- 
tra is not high enough for abundance determinations. We limit 
ourselves to qualitative statements as the absolute abundances in 
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our models are sensitive to uncertainties in the initial composi- 
tion and the chemical profile inside the star. 

6.2.1 . Comments on individual fits: x 2 min < 6.5 

The seventeen systems for which we found fits that we deem ac- 
ceptable are discussed individually below. We focus the discus- 
sion on the characteristics that distinguish the best-fitting mass 
transfer efficiencies. 

26 This is one of the brightest systems, with an 09.5 primary. 
|H03| mention the high quality of the spectra. The asymmetric 
light curve shows evidence for an accretion stream. 
Our non-conservative models (/3 < 0.5) give better fits to this 
system. In fits to conservative models, the primary mass and 
temperature prefer tracks with small initial masses, but they 
cannot be fitted at the same time as the mass ratio and accre- 
tor radius, which are better explained by models with higher 
initial masses. The substantial value of v 2 . is dominated by 

^ mm J 

the contribution of the accretor radius, which is somewhat 
under-predicted in all models. 

This is the semi-detached system with the widest orbit, both 
currently and initially. As a consequence it experienced mass 
transfer in a late stage of the primary's main-sequence phase, 
and the current slow mass-transfer phase will be relatively 
short. 

28 This system, consisting of two O stars, is the most massive 
in the whole sample. It fits well to models with /3 = 0.5 - 1.0 
and an initial mass ratio of about 1.3, although the effec- 
tive temperature is overestimated by more than one sigma 
assuming conservative evolution. According to these mod- 
els the system is currently seen in the early stages of the 
slow mass-transfer phase. For (3 = 0.25 the observed pri- 
mary mass (which requires a model with an initially massive 
primary) cannot be fitted well together with the primary ra- 
dius and the mass ratio, which prefer models with a lower 
initial primary mass and more extreme mass ratio. 

3 1 The best fit to this systems is for (3 = 0.5 with a large initial 
mass ratio (1.7). This model narrowly missed contact dur- 
ing rapid mass transfer, and predicts the current system to be 
near the end of the slow mass-transfer phase, not long be- 
fore core H-exhaustion of the initial primary. Conservative 
(/? > 0.75) models also require a large initial mass ratio to 
fit the currently observed one, but fail to account for the 
observed temperature ratio. These conservative models un- 
dergo temporary contact during rapid mass transfer, and it is 
not excluded that this system evolved fairly conservatively 
through a deeper temporary contact phase, which we cannot 
model properly. 

As a result of the advanced evolution stage, our models pre- 
dict that the C and N abundances have reached equilibrium 
at the surface but oxygen has not, and the O abundance de- 
pends quite strongly on the adopted value of /3. 

32 This system has a fairly wide orbit, and the components have 
very similar masses. The only acceptable fit for this binary 
is for highly non-conservative evolution (J3 = 0.25). The fit 
is not very good C^ min = 5.7) as the models over-predict 
the mass ratio and under-predict the primary mass, but these 
problems are much worse for higher mass-transfer efficien- 
cies. It is possible that a value of f3 < 0.25 would fit this 
system better. 

33 This system fits best for 0.5 < f3 < 0.75, although the best fit- 
ting models overestimate the effective temperature by more 
than 1 sigma. Conservative models overestimate the temper- 



atures even more, whereas models with f3 = 0.25 under- 
predict the masses, but these models are not strongly ex- 
cluded. 

For decreasing values of fi the best-fitting models start with a 
more massive primary and a more extreme mass ratio, so that 
the initial secondary mass is always « 10 M Q . The less con- 
servative models fit this system in a later phase of evolution, 
when more mass has been transferred. Abundance measure- 
ments of C and N could possibly distinguish between a con- 
servative and non-conservative solution, since the predicted 
model abundances for this system depend quite strongly on 

36 This system has a fairly large mass ratio. The best-fitting 
models assume conservative mass transfer /3 > 0.75. 
However, the model temperatures are larger than observed 
while the masses are smaller (a problem that also occurs in 
many of the poorly-fitting systems, Sect. 16.2.2b . The main 
problem with the non-conservative models is that the masses 
are under-predicted even more by the models, by about 2M 

37 With an orbital period of 4 days, this system is one of the 
wider systems in the sample. Only very non-conservative 
models (J3 < 0.5) can fit both the large radius of the Roche- 
lobe filling star and the much smaller radius of its companion 
at the same time as the mass ratio of the system. Like sys- 
tem 26, this binary underwent mass transfer fairly late on the 
main sequence and the primary is currently close to hydro- 
gen exhaustion in its core. 

38 This system fits well to non-conservative models (J3 < 0.53) 
not long before the donor will leave the main sequence. For 
the more conservative models, the observed masses and tem- 
peratures would fit better against initially less massive mod- 
els, while the radii and temperature ratio would fit better 
against more massive models. Also for this system the best- 
fitting models suggest that abundance measurements of C 
and N (and perhaps O) could possibly constrain f3 further. 

41 A good fit to this system, which has a large mass ratio and 
temperature ratio, is obtained for fairly conservative evolu- 
tion (J3 > 0.75), assuming that the system started with a 
mass ratio near unity. The non-conservative models cannot 
fit the observed masses (which prefer models with larger ini- 
tial masses) at the same time as the temperature ratio and the 
radius of the Roche-lobe filling star. 

42 This system can be fitted well assuming non-conservative 
evolution, particularly for j5 — 0.5. The main problem with 
fitting a conservative model is the large radius of mass- 
accreting star (nearly equal to that of its companion) in com- 
parison with its mass and temperature. This system is very 
similar to system 38 in terms of its inferred initial mass and 
mass ratio, but started with a smaller period. 

44 This system is fitted well by our models without a clear pref- 
erence for a certain mass transfer efficiency. It has one of the 
shortest periods in the sample, and in all models mass trans- 
fer started soon after the onset of hydrogen burning and will 
eventually lead to contact (case AS evolution). The smaller 
the value of /3, the closer to contact the system currently is, 
according to our models. 

46 This is the least massive of the semi-detached systems and 
it has the most extreme mass ratio (3.1) of the sample. The 
non-conservative models cannot explain this extreme mass 
ratio, to which we obtain the best fit with a completely con- 
servative model with a large initial mass ratio. This model 
overestimates the temperature ratio, however, and this is the 
major contribution to the rather poor = 6.5. The primary 
temperature is also not perfectly reproduced but the overall 
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fit is acceptable. This model went through a brief contact 
phase during rapid mass transfer. The Roche-lobe filling star 
is fitted right at the end of the slow mass transfer phase, and 
also the accreting star is very close to the end of the main se- 
quence according to this model. The system is on the verge 
of detaching, before entering another phase of case AB mass 
transfer. 

47 This is one of the more massive and hottest systems in the 
sample, and the light-curve solution indicates it is a near- 
contact binary in which the primary star fills 99 per cent of 
its Roche lobe. We find a fairly good fit to a non-conservative 
(J3 = 0.25) model in contact, which however somewhat over- 
estimates the masses and underestimates the temperatures. 
This solution has a very large initial mass of 28 M Q . For the 
more conservative models the observed masses, which prefer 
initially even more massive models, cannot be fitted simulta- 
neously with the other observables. 

48 This system can be fitted well assuming 0.5 < /3 < 0.75. 
More conservative models overestimate the observed small 
mass ratio. The system is currently seen shortly after the end 
of the rapid mass transfer phase. 

49 This is the shortest-period binary in the sample, for which 
the light-curve solution indicates a deep-contact system. A 
good fit solution is obtained for all values of /3; only the most 
non-conservative grid models cannot explain the tempera- 
ture ratio well. In all best-fitting models the system is still 
semi-detached but the accreting star is very close to filling 
its Roche lobe. However, contact is not reached until sub- 
stantially more mass is transferred. 

50 The system can be fitted well assuming non-conservative 
evolution (/3 < 0.5) and an initial mass ratio near unity. If 
conservative mass transfer is assumed the primary mass and 
temperature cannot be fitted well at the same time as the very 
large radius of the mass-accreting star. On the basis of the 
best-fitting model we interpret this system as being in the 
last stage of the slow mass-transfer phase. It is very simi- 
lar to system 42 in terms of its initial parameters and current 
evolution state, except for a smaller initial primary mass (and 
consequently larger age). 

5 1 The alternative detached light-curve solution of this system 
does not fit well to the models, as explained in Sect. 16.1.11 
Conservative models with f3 > 0.75 give good fits to the 
semi-detached configuration, although the masses are un- 
derestimated and the radius of the accreting star is overes- 
timated. This is worse for non-conservative models, which 
also overestimate the radius of the accreting star. The best 
fit is found for a model with a large initial mass ratio, which 
went through a contact phase during rapid mass transfer and 
is currently undergoing case AB mass transfer. The C/N ratio 
of the best-fitting models decreases with so that its deter- 
mination might further constrain the mass transfer efficiency. 

6.2.2. Comments on individual fits: x 2 ^ m > 10 

For the systems, for which we find no acceptable fits, we indicate 
which of the three common conflicts occur (see the start of this 
subsection). We divide the systems in two groups, based on their 
X 2 values. We start with the eight systems that have 10 < x 2 ■ < 
20: 

22 This system shows the ^-conflict. Only the conservative 
models can match the very large temperature ratio of this 
system, as well as the moderately small value of R^/Ra- The 
mass ratio is then over-predicted by the models. 



23 This system exhibits the q- and Mr-conflict. The best fitting 
model, obtained assuming — 0.5, fits the observed radii but 
underestimates the mass of the donor star by about 3 M . 

34 The /^-conflict and a moderate form of the g-conflict occur 
for this system. The best fit is obtained with conservative 
models. 

39 The q- and /^-conflicts are present as well as a moderate in- 
verse Mr-conflict. 

40 A combination of all three conflicts occur for this system. 

43 This system, which shows the ^-conflict and a moderate in- 
verse Mr-conflict, has one of the highest effective tempera- 
tures. In the best fit, assuming /3 = 0.5, the temperature ratio 
is under-predicted and the mass ratio is over-predicted by the 
models. 

45 For this system with a ^-conflict, the conservative models 
provide a reasonably good fit to nearly all observed parame- 
ters for this system except the mass ratio. 

52 The alternative detached solution for this model could not 
be fitted well with our models. The same is the case for 
the semi-detached solution. The system shows the q- and in- 
verse Mr-conflict. The radius of the accretor is also under- 
predicted. 

Six systems have very poor fit solutions, with^ nin > 25: 

24 This systems shows a severe form of the ^-conflict and can- 
not be matched for any assumption of p. 

25 This system suffers from a severe g-conflict as well as the 
^-conflict and the inverse Mr-conflict. 

27 This system shows the q- and /^-conflict. 

29 This system exhibits the q- and Mr-conflict. 

30 This system shows a severe form of the ^-conflict. The large 
X 1 • is dominated by the contribution of the extreme tem- 

^ mm J 

perature ratio, which is the largest in the sample. The other 
two common discrepancies, the R- and Mr-conflict, are also 
very serious in this system, Rd/Ra is even slightly below 
unity. The photometric flux ratio is incorrectly reported in 
H05 as 1.2; it should be 2.7. The spectral type of the sec- 
ondary component is B2 or later. 

35 This system is the least massive of all semi-detached systems 
and cannot be fitted well because of the ^-conflict. 

6.3. Summary of fit results 

As we have argued in Sect. [5] due to the nature of our fit method 
the absolute value of x 2 oin is not a very good measure of the 
goodness-of-fit. The observables we have fitted are not truly in- 
dependent, and their errors may not in all cases be representative 
of the true uncertainties. This applies especially to the tempera- 
ture ratio, which is determined from a flux ratio that is not very 
sensitive to temperature for OB stars. Also the error on the pri- 
mary temperature is only a rough estimate. Therefore we feel we 
can be fairly generous in accepting fit solutions with large ;kfjL in , 
and we have argued in the previous subsections why we regard 
fits with^ in < 6.5 as acceptable. 

In Fig. [6] we show the cumulative distribution of^ in for the 
detached and the semi-detached systems. For about 80% of the 
detached systems and only about 55% of the semi-detached sys- 
tems we obtain fits with^ in < 6.5. Interestingly, both distribu- 
tions show a break between this and larger values of x 2 - There 
is no reason to expect that, if any of the observational errors 
were underestimated, this should occur preferentially for semi- 
detached binaries. It is noteworthy that the ^-conflict, which af- 
flicts a large fraction of the semi-detached fits, only occurs in 
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Fig. 6. Fraction of detached (full line) and semi detached sys- 
tems (dotted line) for which a fit solutions were obtained with 
^-^ in lower than a certain value. We have excluded system 1 from 
the detached systems and system 52 from the semi-detached sys- 
tems. 
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Fig. 7. The best fitting mass transfer efficiency ft is plotted 
against the best fitting initial orbital period P, . The symbols indi- 
cate the magnitude of^j _. The error bars indicate the l<x range 
of the two fit parameters; if the best-fitting p = 0.25 the error bar 
is extended downward by an arbitrary amount. No strong trend 
is found but there is a hint that systems with small initial orbital 
periods evolve more conservatively than wider systems. 

one detached system (11). We must therefore conclude that our 
models do significantly worse at explaining post-mass transfer 
binaries than detached systems which have not yet undergone 
mass transfer. We will discuss possible causes and the limita- 
tions of our models in Sect. [8] 



Table 4. For each semi-detached system x^. is given as func- 
tion of p. For each system the value corresponding to the best fit- 
ting p is marked. The last column indicates the estimated range 
of p which fits well for a certain system. 
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Fig. 8. As Fig. |7]y6 against mass ratio q,. 



7. Conclusions on the mass transfer efficiency 

Although a large fraction of the semi-detached systems cannot 
be fitted well with any assumption for the mass transfer effi- 
ciency, we will focus on those systems for which we do obtain 
reasonable fits to see if we can learn anything about the effi- 
ciency of mass transfer from the comparison of models and ob- 
servations. In Table R]y 2 . is given as a function of the mass 
transfer efficiency. The last column gives an estimate for the 



range of p for which a certain system can be fitted, as explained 
in Sect. [5] Some of the semi-detached systems fit significantly 
better against conservative models (p > 0.7), for example sys- 
tem 36, 41, 46 and 51, while other systems show a strong pref- 
erence for the non-conservative models (J3 < 0.5), for example 
system 32, 37 and 47. There is no single value of p for which 
we can fit all systems. From this we conclude that the simplified 
assumption of conservative evolution is not valid. Neither is the 
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Fig. 9. As Fig. [7]for /3 against primary mass Mij. 

often made assumption of/3 - 0.5. Binary systems with different 
initial parameters cannot be described with one fixed value of /?. 
The spread in best-fitting /3 values covers the full range of /3. 

In order to understand the spread in /3, we searched for cor- 
relations between the preferred /3 and the initial parameters. In 
Figure [7] the preferred /3 is plotted against the initial orbital pe- 
riod. The correlation is not strong but there is a region in the 
upper right part of the diagram where we find no systems: the 
initially wider systems prefer less conservative models. If real, 
this correlation might be explained by the interplay of two mech- 
anisms: spin up of the accreting star and tidal interaction. The 
mass-accreting star is spun up as it accretes angular momentum. 
Tidal interaction tends to synchronize the stars with the orbit 
and prevent the accreting star from reaching critical rotation. In 
wider systems, where tidal interaction is less efficient, the ac- 
creting star rotates faster and may lose mass along its equator. 

One can expect a strong correlation of the mass transfer ef- 
ficiency with initial mass ratio, because in systems with extreme 
mass ratios the accretor will expand rapidly such that a con- 
tact system is formed. Mass and angular momentum can than 
be lost from the outer Lagrangian point. However, we found 
no correlation of the mass transfer efficiency with the mass ra- 
tio in the expected direction, see Figure [8] This might be be- 
cause the range in current mass ratios available in the sample 
is small, as the sample is strongly biased towards systems with 
equal masses. No significant correlation with the initial primary 
mass was found either, see Figure[9] 

8. Discussion 

It is striking that we can fit the detached systems in the sample 
quite well, while a large fraction of the semi-detached systems 
cannot be explained with our models. We seem to understand 
the evolution of detached systems well, i.e. essentially the evo- 
lution of two single stars, while our models fail to account for a 
large number of systems in slow mass transfer, especially their 
temperatures. 

In our models we assume a constant mass transfer efficiency 
throughout the whole evolution. However, mass loss from the 
system is more likely to occur mainly during the rapid phase of 
mass transfer. A more realistic description of the mass transfer 
phases is to be preferred. For instance, spin-up by accretion and 
associated rotatio n-induced mass lo ss ca n be taken into a ccoun t, 
as was done by IWellsteinl d2001l) and iPetrovic et ail d2005l) . 



However, this treatment introduces other uncertainties, e.g. for 
the strength of rotation-induced mass loss and the timescale for 
tidal interaction, which are both uncertain. 

For non-conservative evolution we assume that the lost mass 
carries the specific orbital angular momentum of the accreting 
star. While the underlying assumption (that the accreting star 
ejects the surplus of mass isotropically) is reasonable, it may 
be worthwhile to test the effect of different assumptions. For in- 
stance, somewhat stronger angular momentum loss than we have 
assumed here may improve some of the fits exhibiting the R- 
confiict. 

Our models with extreme initial mass ratios show that some 
systems undergo a temporary contact phase, after which the sys- 
tem restores thermal equilibrium and becomes semi-detached 
again. A significant fraction of our models with extreme initial 
mass ratios fail to converge during rapid mass transfer and we 
stop our calculations if deep contact occurs, as our code contains 
only a crude model for these contact situations. Improvement of 
the description of the contact phase might enable us to fit more 
systems properly. 

We assume the Schwarzschild criterion for stability against 
convection. As a result mass-accreting stars always rejuvenate 
(i.e. increase the size of their convective core). The question of 
whether stars rejuvenate or not depends on the still poorly under- 
stood efficiency of semi-convection. It is possible that the steep 
molecular weight gradient, built up before th e onset of mass 
transf er, prevents the star from rejuvenating dBraun & Langerl 
fl995l) . This shortens the remaining lifetime of the accreting star 
and affects its position in the H-R diagram as a function of time. 
Whether this effect would improve any of the fits remains to be 
investigated. 

In principle mass transfer can lead to extra mixing, for exam- 
ple because of spin up and resulting differential rotation, which 
is not included in our models. If the extra mixing is significant, 
then chemical abundance profiles in the star will become flatter. 
In general this results in more compact and therefore hotter stars. 
As this only affects the mass-accreting star, which is the hotter 
component, it will lead to more extreme temperature ratios. This 
might help to solve the discrepancy between the observed and 
modeled temperature ratios. 

To account for the discreteness of our model grid we esti- 
mated the spacing in our grid by the initial spacing. The actual 
grid spacing changes during the evolution and it might even de- 
pend on /3. Interpolating between the models would be a solution 
but this is not simple considering the many dimensions in the 
grid. We investigated the effect of different values for o-± and 
we found no evidence for systematically under- or overestimat- 
ing the actual grid spacing for evolved systems. 

Our fit method might also be improved by using more di- 
rectly measured observables. We recall that the most commonly 
occurring conflicts between our models and the observations are 
related to the temperature ratio and the radius ratio. For this sam- 
ple the temperature ratios were derived from the /-band flux ra- 
tios as determined from the light-curve solution. Rather than fit- 
ting to the temperature ratio, we can use the /-band flux ratio and 
its error directly as an observable. The light-curve solution often 
yields the sum of the radii to a higher accuracy than the ratio of 
the radii, as a result of which the errors in the individual radii are 
correlated. It may be preferable to use the sum and the ratio of 
the radii as observables, rather than the individual radii. 
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9. Summary and outlook 

We computed an extensive grid of detailed binary evolution 
tracks at a metallicity Z = 0.004 for a variety of mass transfer 
efficiencies. Our models are available to the astronomical com- 
munity. We plan to extend the grid to a larger range of binary 
parameters and to other metallicities in the future. 

We have compared our models to observations of 50 double- 
lined eclipsing binaries, by fitting evolution tracks to each in- 
dividual system. For the detached systems in the sample we 
find generally good agreement between observations and mod- 
els. Our models can also explain a large fraction of the semi- 
detached systems although the overall agreement is less good. 
We identified three common conflicts between the models and 
the observations. 

(I) Several semi-detached systems in which the stars have 
fairly equal masses are observed to have a large temperature ra- 
tio, more extreme (in some cases much more extreme) than our 
models predict. (II) In some systems the ratio of donor radius 
to accretor radius, Rd/Ra, is smaller than shown by our mod- 
els, in which case the radius of the accreting star is often under- 
predicted. (Ill) A third problem is the combination of a large 
observed primary mass with a relatively low temperature which 
cannot be matched simultaneously by our models. The inverse 
(small mass together with high temperature) also occurs in a few 
systems. 

Conservative mass transfer cannot explain these case A mass 
transfer systems. We find a large spread in the best fitting mass 
transfer efficiency. There is no single assumption for/? which can 
explain all semi-detached systems. We conclude that the often 
made assumption of a single constant /3 for all case A binaries 
is not valid. We find a hint that initially wider systems fit better 
to non-conservative models. This might indicate the importance 
of spin up and tidal interaction in determining the efficiency of 
mass transfer. 

We have shown that, assuming non-conservative mass trans- 
fer, more systems will avoid contact during rapid mass transfer 
and that the duration of the slow mass transfer phase becomes 
longer. This can in principle give a strong test on the mass trans- 
fer efficiency by comparing the number of detached and semi- 
detached systems. Good understanding of the selection effects 
and the initial mass ratio distribution are needed for such a test. 

Although the observed sample we used is currently the 
largest single set of stellar parameters determined for massive bi- 
naries in any galaxy, it is limited in primary masses, mass ratios 
and orbital periods. It would be interesting to extend this study to 
all currently known semi-detached double-lined eclipsing bina- 
ries in our galaxy and the Magellanic Clouds. This would enable 
us to study a larger range of primary masses. Studying a larger 
range of orbital periods and mass ratios will however be difficult 
as all such samples are biased towards small periods and equal 
mass ratios. 

The largest discrepancies between models and observations 
are seen in the temperature ratios. The /-band flux ratios from 
which these are determined are not very sensitive to temperature 
differences as they cover only the Raleigh-Jeans tail of the spec- 
trum. Light curves in the U and B-band would lead to more reli- 
able determinations of the temperature ratio. From high quality 
spectra accurate temperatures can be obtained by fitting model 
atmospheres. 

Our results indicate that surface abundance determinations 
of nitrogen and carbon can potentially constrain the mass trans- 
fer efficiency strongly. If mass transfer is a non-conservative pro- 
cess, then in order to reach a given mass ratio the donor must lose 



a larger fraction of its mass than for conservative mass transfer. 
Therefore deeper layers of the star are revealed, in which CN- 
cycling has increased the N/C ratio. In Section [6721 we have in- 
dicated several particularly interesting and promising targets for 
further study. In order to determine accurate temperatures and 
surface abundances, high-resolution spectroscopic follow-up ob- 
servations are needed. 
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Fig. 10. See Fig|4]for description of the symbols and the units. The observed parameters are plotted together with the best fitting 
model, see table [2] and [3] For the detached systems we plotted conservative tracks. 
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